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Abstract. We implemented a joint weak measurement of the trajectories of two 
photons in a photonic version of Hardy's experiment. The joint weak measurement 
has been performed via an entangled meter state in polarization degrees of freedom 
of the two photons. Unlike Hardy's original argument in which the contradiction is 
inferred by retrodiction, our experiment reveals its paradoxical nature as preposterous 
values actually read out from the meter. Such a direct observation of a paradox gives 
us new insights into the spooky action of quantum mechanics. 
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Although it is natural to ask what is the value of a physical quantity in the middle 
of a time evolution, it is difficult to answer such a question in quantum mechanics, 
especially when post-selection is involved. Hardy's thought experiment [1] is a typical 
example in which we encounter such a difficulty. Figure [T](a) shows a photonic version 
of the experiment, which was demonstrated by Irvine et al. recently [2]. The scheme 
consists of two Mach-Zehnder interferometers MZl and MZ2 with their inner arms 
(Oi, O2) overlapping on each other at the 50:50 beam sphtter BS3. If photons 1 and 
2 simultaneously arrive at BS3, due to a two-photon interference effect, they always 
emerge at the same port. This corresponds to the positron-electron annihilation in the 
original thought experiment pj. The path lengths of MZl are adjusted so that photon 1 
should never appear at Ci by destructive interference, when photon 2 passes the outer 
arm NO2 and thus has no disturbance on MZl. The path lengths of MZ2 are adjusted 
similarly. Then, a coincidence detection at Ci and C2 gives a paradoxical statement on 
which paths the detected photons have taken. The detection at Ci (C2) implies that 
MZl (MZ2) has been disturbed by photon 2 (1) travelling along O2 (Oi). We may thus 
infer that the conditional probabilities satisfy 



On the other hand, if both photons had taken the inner arms, the coincidence detection 
would have never happened due to the two photon interference. Hence we infer that 



The two inferred statements are apparently contradictory to each other, which is the 
well-known Hardy's paradox. 

One may argue that we should abandon the attempt to address the question itself, 
on the ground that the trajectory of photons cannot be measured without utterly 
changing the time evolution. But this reasoning is not necessarily true if we are 
allowed to repeat the same experiment many times. Aharonov et al. has proposed 
weak measurement O H], in which a measurement apparatus (meter) interacts with 
the system to be measured so weakly that the state of the system is not significantly 
disturbed. The readout of the meter from a single run of experiment may be subtle 
and noisy, but by taking the average over many runs we can correctly estimate the 
expectation value of the measured observable, {ip \ A\ ip), when the initial state of the 
measured system is In this setup, we may ask what is the averaged readout over 
the runs in which the system is finally found to be in a state | 0). In the limit of no 
disturbance, this gives an operational way of defining what is the value of A at the 
middle of a time evolution from \ ip) to \ (p) , and is found to be given by the real part of 
the following expression 



O El [71 [HI [9l [TOj [TT] and experimental observations of weak values were reported 
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Figure 1. Experimental setup for Hardy's paradox, (a) The photonic version of 
Hardy's thought experiment. Each of the MZ interferometers MZl and MZ2 is 
composed of four 50:50 beam sphtters (BS). (b) The schematics of the experiment 
with joint weak measurement. Entangled photon pairs are generated via spontaneous 
parametric down-conversion from a pair of type I phase matched 2mm thick BBO 
crystal T7J pumped by a UV pulse (a central wave length 395nm and an average 
power 180mW). The UV pulse is taken from the frequency doubled Ti:sapphire laser 
(wavelength, 790nm; pulse width, 140fs; repetition rate, 76MHz). Quartz crystals are 
used to compensate the group velocity mismatch and adjust the relative phase between 
horizontal and vertical polarization state. Polarization dependent phase shifts in MZs 
are compensated by quartz plates in the MZs so that MZs give the same phase shifts 
to any polarization state. Extra laser diodes (LDs) and photo diodes (PDs) are used 
for adjusting and stabilizing the optical lengths of MZs via piezo stage (Ps). The 
photons are detected by photon detectors (Cs). The observed visibility of two photon 
interference at BS3 was 97.8 ± 0.3% for horizontally polarized photons. 



Suppose that weak measurements of trajectories are applied to Hardy's experiment 
at the shaded regions in Fig. [T](a). The state of the photons entering these regions 
is = {\N0i)\02) + |0i)|A^02) + \NOi)\N02))/V3, and the coincidence detection 
retrodicts the state leaving the regions to be |0) = (|A^Oi) - |0i))(|A^02) - |02))/2 [S]. 
Then the weak values can be calculated to be 

|Oi,02)(Oi,02|w = 0, \NO,,N02){NOuN02U = -1, 
|Oi,iV02)(Oi,iV02|w = 1, |iVOi,02)(iVOi,02|w = 1. (4) 

The first equation implies Eq. ([2]) holds. We also see that Eq. ([T]) holds since, for 
instance, \0i){0i\^ = |0i, 02)(0i, 02|w + |C'i, A^02)(0i, A^02|w = 1- Hence the readout 
of the meter is indeed consistent with both of the naively inferred conditions ([H) and ([2]). 
The reason why these two contradictory conditions are satisfied at the same time can 
now be ascribed to the appearance of a negative value, |A^Oi, N02){N0i, A^02|w = — 1- 
It implies that the average readout over post-selected events falls on a value that never 
appears if no post-selection is involved. 

In this paper, we report an experimental demonstration of weak measurements 
applied to Hardy's experiment. Since the observables to be measured are of path 



correlations between the two photons, we need to implement a joint weak measurement 
that can extract such correlations. Although a number of schemes for joint weak 
measurement have been proposed [I8l|19], here we propose the most intuitive approach 
using meter qubits initially prepared in an entangled state. We conducted measurements 
with varied strength, and confirmed that it works properly when each photon is injected 
to a fixed arm. We also measured the visibility of the interferometers to verify that the 
disturbance vanishes in the limit of weak measurement. Then we conducted Hardy's 
experiments and observed that, as the measurement became weaker, three readouts 
moved toward the conditions ([1]) and ([2]) while another readout went to negative. 

Let us first introduce our scheme of joint weak measurement for two qubits in a 
general setting. Let {| 0), 1 1)} be the standard basis of a qubit. For weak measurement 
of the observable | 0)(0 | of a qubit, one can prepare another qubit in a suitable state 
and apply a controUed-NOT (C-NOT) gate between the two qubits [H]. Here we want 
to carry out weak measurements of observables \ kl)sj^s2{kl \ {k,l = 0,1) for the two 
signal qubits, si and S2, and we may use two meter qubits, mi and m2, and apply two 
C-NOT gates in parallel. The signal (meter) qubit corresponds to the control (target) 
qubit and the meter qubit is flipped when the signal qubit is in |1). Note that in our 
implementation of Hardy's experiment, the photon path corresponds to the signal qubit 
and its polarization to the meter qubit. We found that the desired weak measurements 
are achieved by a simple choice of the initial state of the meter qubits, 

+ e(|01) + |io) mim2 

= {6 - e)\m)rmm2+ £^\kl)rmm2: (5) 

kl 

where 5^ + ?>e'^ = 1 and 6 > e > 0. When the signal qubits are in \ip)s = Cjj|ij%^s2, 
the application of the parallel C-NOT operation Ugm results in 

Usm\^)s\Om 

~ ^ ] ^ijl'^j) siS2[{^ ^)l^j)mim2 ~^ ^ ^ ^ |^0m.ir?i2] 
ij kl 

= ^[((5 - e)\kl)s,s2{kl\ + e]\ij)s\kl)rr,,m2- (6) 

kl 

The meter qubits will then be measured in the basis {| kl)mim2} to produce an outcome 
{k, I). If the signal qubits are in state \ij)siS2 initially, the probability Pm{k, I) of outcome 
{k,l) becomes 

'"^'^"l^^ ((A:,0^(.,J))- 

We see that the outcome (i, j), corresponding to the state of the signal, has a larger 
probability than the others. The contrast 5^ — e"^ will be regarded as the measurement 
strength. When 5^ — = 1 {6 = 1, e = 0), this scheme gives a projection measurement 
(strong measurement) on the signal qubits. Decreasing 6"^ — makes the measurement 
weaker, and at 5^ — = (5 = e = 1/2), the operation of Ugm introduces no disturbance 



on the signal qubits (Appendix A). Let us introduce a normalized readout R{k,l) by 

R{k,l)^[P^{k,l)~e']/{5'-e'). (8) 

Then, regardless of the value of 5"^ — > 0, we have R{k,l) = 1 for {k,l) = and 
R{k, I) = otherwise. When the initial state of the signal qubits is | 4')s = %Kj)sis2) 
the normalized readout becomes R{k,l) = |cfc^/p, which coincides with the expectation 
value (I kl)s-^s2{k^ I)- Of course, for the estimation of R{k,l) and hence of Pm{k,l), we 
need to repeat the preparation of the initial state | and the measurement many 
times. 

Now in a situation like Hardy's experiment, we are interested in the readout of the 
meter on condition that the signal qubits are finally measured to be in state \(f))s- From 
the definition (JSj), we have 

s{<P\Usm\lp)s\Om 

= 5^[(5-£)|A;O.,s.(fc/|w + ^].(0|^).|A;On^. (9) 

kl 

The conditional probability of the outcome {k, I) is then given by 

PUk, m = 1(5 - e)\kl)s,s,{klU + - {S - ef(] (10) 

with ^ = 1 — I lu)siS2(^ilwP- The normalized readout is found to be 

R(k,l\4>)^[PUkJ\4>)-e']/{6'-e') (11) 
^ 2e{d - e)Re[\kl)s,sAklU] + {S - £)^[| |A;/),,,,(fc/|„P + e^C] 

{6^ - e^)[l - {6 - en] 
^Re[\kl)s,s,{klU] {6'-e^^0). (12) 

We see that the normalized readout in the limit of no disturbance coincides with the 
real part of the weak value. 

In our experiment, the signal qubit corresponds to the paths taken by the photon 
as I 0) = I NO) and | 1) = \0), whereas the meter qubit is assigned to the polarization 
of the same photon as I 0) = I +) and 1 1) = | -), where |±) = {\H) ± \V))/V2 and H/V 
is the horizontal/vertical polarization. The interferometers in Fig. [T](a) are constructed 
so that there is no polarization dependence, except for the shaded region where wave 
plates are placed to reahze Usm- The polarization dependence of the coincidence events 
is then analyzed to determine Pm{k, ^|0). 

The detail of the setup is shown in Fig. [TJ^b). We generate the photon pairs in a 
nonmaximally entangled state ri\HH)i2 + f]\VV)i2 from PDC, where 7f + ff = 1 and 
7] and f] are properly adjusted to be real numbers. The coefficients rj and f/, which 
can be tuned by HWP3, are determined by the ratio between the coincidence counts 
of \HH) and |^^). The state is transformed into \C)m by rotating the polarization 
of two photons by half wave plates HWPl and HWP2. The measurement strength is 
simply calculated by the measured r] and f) and the angle of HWPl and HWP2. After 
the photons pass through BSl, BS2 and BS3, the two-photon state is represented as 
\'^)s\C)m- The C-NOT operations between signals and meters are simply performed by 
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Figure 2. Relation between the visibility of MZl (2) and the measurement strength. 
The solid curve represents the theoretically expected visibility, 2e{S + e). 



flipping the polarization as |±) |=f) via HWP4 and HWP5 in the overlapped arms, 
Oi and O2. The polarization basis (| +), | — )) for the meter measurement is selected 
by HWP and PBS just before detectors Ci and C2. The coincidence counts, typically 
a few thousands, are accumulated over 5-12 min per basis to determine each R{k,l\(j)). 
All error bars assume the Poisson statistics of the counts. 

First we show the observed trade-off between the measurement strength and the 
visibilities of interference at MZl and MZ2. The measurements were done as follows: 
We blocked the arm 02(0i) and recorded the coincidence detection at Ci and C2 by 
changing the path length via piezo stage Pi(P2). In order to measure the polarization- 
independent visibilities, we measured the interference fringes in the states | + +), | H — ), 

I — h), and I ) independently and calculated the average visibility. The visibilities for 

various measurement strengths are shown in Fig. [2l We can clearly see that weakening 
the measurement makes the visibility larger (Appendix B). 

We also show that our measurement gives proper values of the readouts R{k,l\(j)) 
when the photons travel the fixed arms \Oi,N02),\NOi,02), and |A^0i,A^02). We 
blocked the arms of no interest and determined R{k, from the coincidence detection 

on the basis { | + +), | H — ), | — h), | ) } . The experimental results are plotted 

in Fig. [3l In any measurement strength, the corresponding readout for the signal state 
is close to 1 and the other readouts are close to 0. This clearly shows the apparatus 
properly measures the trajectories of the photons. The systematic deviation from the 
expected values at small values of 5^ — is believed to be the results of very tiny 
misalignments of the axes of the wave plates. Such an enhancement is a feature of weak 
measurement and can be used for super- sensitive measurements as shown in [T6] . 

With this meter, we performed measurements of the trajectories of the photons in 
Hardy's thought experiment. The observed readouts are shown in Fig. |H Let us see 
the results around 0.3 measurement strength. The visibility above 0.9 shown in Fig. [2] 
suggests that the measurement is fairly weak, and the systematic deviation in Fig. [3] is 
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Figure 3. Experimental results when the initial states are in {a,)\Oi, NO2), 
{h)\NOi,02), and ic)\NOi, NO2) . 
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Figure 4. Observed readouts in Hardy's experiment. The curves that roughly fit the 
data points are of a simple model of imperfections, and arc only shown to guide the 
eye. 



small at this strength. For each of the photons, the readouts suggest a high probability 
of taking the overlapped arm: R{Oi\(f)) = R{Oi, N02\(j)) + R{Oi,02\(p) = 0.72±0.09and 
i?(O2|0) = 0.53±0.09. Nonetheless, the readout for the joint probability is around zero, 
namely, R{Oi, O2I0) = -0.18 ± 0.05. We also see that R{NOi, N02\(p) = -0.43 ± 0.04 
takes a large negative value, which cannot be interpreted as a probability. Although 
the statistical and systematic errors are large around zero measurement strength, we 
see that i?(Oi, A^Oal^), i?(iVOi, O2I0), i?(Oi,O2|0), and R{NOi, N02\(l)) approach the 
expected weak values in Eq. (jl]) when the measurement strength goes to 0. 

Let us emphasize here that our experimental results by themselves can elucidate 
the paradoxical nature of Hardy's experiment, without any reference to the theory 
explaining how our measurement works. From the results in Fig. [3], we see that our 
measurement apparatus should properly give the probabilities on the trajectories of the 
photons, if such quantities ever exist. Then, in Hardy's setup, the same measurement 
presents a contradictory statement in Fig. H] that the probability of finding a photon 
in arm Oi approaches 1 and so does the probability for arm O2, whereas the joint 



probability of photons being in both of the arms stays at about 0. Moreover, the readout 
points to a negative value for the joint probability for arms NOi and N02- Unlike 
Hardy's original argument, our demonstration reveals the paradox by observation, rather 
than inference. 

We have experimentally demonstrated a joint weak measurement with an entangled 
photon pair and directly observed paradoxical results in Hardy's thought experiment. 
Our demonstration clearly reproduces Hardy's paradox in an experimentally accessible 
manner. Weak measurements have attracted attention due to the possibility of achieving 
small-backaction and high-sensitivity measurements by simple optical setups. We believe 
the demonstrated joint weak measurement is useful not only for exploiting fundamental 
quantum physics, but also for various applications such as quantum metrology and 
quantum information technology. 

Note added. During the preparation of this manuscript, we came to know of a 
work by Lundeen and Steinberg [20], which has demonstrated similar results by using, 
interestingly, different methods for joint weak measurement and for constructing Hardy's 
interferometer. 
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Appendix A. 

The positive operator valued measure (POVM) for the measurement on the signal qubits 
is given by 

Tiki = 6^\kl),,sAkl\+e^ Yl l^^%i^2(ul 

ij^kl 

= {S^-e'')\kl),,s,{kl\+e'. (A.l) 

The same POVM can be realized also by a nonentangled initial state of the meter qubits, 
5'|00)^,^,(00| +£2(|01)^^^^(01| + |10)^,^,(10| + |11)™(11|). This separable state of 
the meter qubits, however, fails to derive a joint weak value, because the signal qubits 
are utterly disturbed even if 6^ — = 0. 

Appendix B. 

A general relation between the visibility V and the measurement strength was discussed 
in [21]. It was shown that V"^ + < 1 holds, where K = 1 — 2pe„ stands for 
the measurement strength defined in terms of the error probability Perr in the which- 
path measurement. In our case, the initial marginal state of the meter qubit mi is 



a mixed state, which is, from Eq. (j^]), calculated to be a mixture of state | 0o)mi = 

(5| 0)mi + ^1 l)mi)/ v^Po with probability po = + and | 0i)mi = 

with pi = 2e^. The state | 0o)mi leads to a measurement strength Kq = / 

and visibility Vq = ^^/T^^l{^ = 2e6 /{S'^+e'^), whereas | 0i)mi leads to Ki = and Vi = 1. 

Hence, on average, the measurement strength becomes K^ve = PoKq + PiKi = — 

and Vave = PoK) + Pi^i = 2£:(5 + e), which is shown as the solid curve in Fig. 3. Due 

to the mixture, the averaged values Kave and V^ve do not saturate the general bound 
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